Hafnium binary alloys from experiments and first principles 
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Despite the increasing importance of hafnium in numerous technological applications, experimen- 
tal and computational data on its binary alloys is sparse. In particular, data is scant on those binary 
systems believed to be phase separating. We performed a comprehensive study of 44 hafnium binary 
systems with alkali metals, alkaline earths, transition metals and metals, using high-throughput first 
principles calculations. These computations predict novel unsuspected compounds in six binary sys- 
tems previously believed to be phase separating. They also predict a few unreported compounds 
in additional systems and indicate that some reported compounds may actually be unstable at low 
temperatures. We report the results for the following systems: AgHf, AlHf, AuHf, BaHf*, BeHf, 
BiHf, CaHf , CdHf, CoHf, CrHf, CuHf, FeHf, GaHf, HfHg, Hfln, Hflr, HfK*, HfLa*, HfLi*, HfMg, 
HfMn, HfMo, HfNa*, HfNb*, HfNi, HfOs, HfPb, HfPd, HfPt, HfRe, HfRh, HfRu, HfSc, HfSn, HfSr*, 
HfTa*, HfTc, HfTi, HfTl, HfV*, HfW, HfY*, HfZn, and HfZr 

(* = systems in which the ab initio method predicts that no compounds are stable). 



High-throughput (HT) calculations of material proper- 
ties based on density functional theory (DFT) have been 
developed in recent years for theoretically guided mate- 
rial discovery and improvement [l], H, H, m H, H, 0, @|- 
These calculations give insights into trends in alloy prop- 
erties and indicate possible existence of hitherto unob- 
served compounds. In this paper we apply the HT ap- 
proach to a comprehensive screening of hafnium inter- 
metallic binary alloys. This choice is motivated by the 
wide array of technological applications of hafnium alloys 
in contrast with their scant theoretical discussion in the 
literature. 

Hafnium is primarily used in the control and safety 
mechanisms of nuclear reactors, because of its high cross- 
section for neutron absorption and its high corrosion re- 
sistance @. Hafnium cladding of nuclear fuel rods is 
expected to be an important element in the design of fu- 
ture advanced reactors [l(j | . Hafnium is used extensively 
as an alloying element in nickel-, niobium-, and tantalum- 
based superalloys, which are designed to withstand high 
temperatures and pressures. It is an important addi- 
tion to some titanium, tungsten and molybdenum al- 
loys, where it forms second-phase dispersions (with car- 
bon) that improve material strength under extreme con- 
ditions H, [ll|. Hafnium alloys are also used in medi- 
cal implants and devices, due to their biocompatibility 
and corrosion resistance (see for example [Hi])- Nickel- 
titanium-hafnium alloys exhibit shape memory behav- 
ior with high martensitic transformation temperatures 
and good mechanical properties [HI]. Hafnium is added 
to aluminum-magnesium-scandium alloys, widely used in 
aerospace applications, to increase their strength follow- 
ing high temperature thermomechanical processing [14| . 
Some intermetallic compounds of Hf and the transition 



metals Fe, Co, Pd and Pt have been investigated as 
hydrogen-storage materials because of their capability to 
form hydrides with high hydrogen to metal ratios at room 
temperature [T^. 

Hafnium oxide based compounds have recently found 
wide application replacing silicon oxide as high-fc di- 
electrics in the production of integrated circuits (l6| . 
This has motivated a few first principles studies of the 
dielectric properties of hafnium silicates (see for example 
[l7|). Perovskite alkaline metal hafnates (e.g. CaHfC>3, 
BaHf03 and SrHfOs) have been investigated for vari- 
ous optical and electronic applications as well as sub- 
strates for perovskite superconductor films, due to their 
high stability in the fabrication process of such films and 
their small crystallographic mismatch (see e.g. [HI). Al- 
though zirconium has been tried as a cheaper substitute 
for hafnium, in most of these applications only hafnium 
produces the desired properties [9j. 

Despite this wealth of existing and potential applica- 
tions, computational studies of hafnium compounds are 
few. In most of these studies, just a few specific struc- 
tures have been investigated 0, H HH S3, [|l H [H, US 
[H H, [H, [H HH pi pE H . In only three cases, Al-Hf , 
Cu-Hf and Hf-Nb [I9l.l35l.r36l. H3] , have a large number of 
structures been examined. In the Al-Hf system, the en- 
ergies of 22 intermetallic structures have been calculated, 
four of which are ground states and constitute the con- 
vex hull of the system [3^, HH (The convex hull connects 
all of the ground states on a concentration vs. enthalpy 
diagram. The enthalpies of all other structures lie above 
the tie lines between the ground states. Thermodynam- 
ically, the convex hull represents the Gibbs free energy 
of the alloy at zero temperature [H ) . In the Cu-Hf sys- 
tem, energies of 28 structures have been calculated and 
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four identified as ground states [37j . In both cases, differ- 
ences have been found between the phase diagrams based 
on calorimetric measurements and the calculated ground 
states. In contrast, positive formation enthalpies were 
calculated for 13 Hf-Nb intermetallic structures [l9j], in 
agreement with experimental observations that the sys- 
tem is phase-separating. 

In several systems, the stability of a few competing 
structures have been examined but the complete con- 
vex hull of the intermetallic system has not been ad- 
dressed. These include the ground state properties of 
five different structures of Hf 3 Mo, HfMo and HfMos J2I 
and the three Laves phases of HfCr2 [21(, HfFe2 [2al23l|. 
HfV 2 [Mj2land HMn^JpJ. Single specific structures of 
Hf 2 Fe 12013], Hf 2 Co [lf7Hffl2 [M 111, Hf 2 Ni [13, 
HfTe 5 [H, HfRu, HfRh, Hflr and HfPt[34| have also 
been studied. Some of these compounds were studied be- 
cause they are used in hyperfine field measurements by 
the time dependent perturbed angular correlation (TD- 
PAC) technique, for which 183 Hf is an excellent probe, 
that can be easily compared to first principles calcula- 
tions of the local electric field gradient in the vicinity of 
the probe. 

In this paper we report a comprehensive computational 
study of the low temperature stability of 44 hafnium bi- 
nary systems; hafnium combined with alkali metals, the 
alkaline earths, transition metals and metals. We apply 
the ab initio formation energy criterion which has been 
shown to be reliable for binary systems [!, 0, HI] . The 
calculations were performed using the high-throughput 
framework AFLOW [U, 13^ ] , employing ab initio calcula- 
tions of the energies with the VASP software [4cj ]. We 
used projector augmented waves (PAW) pseudopoten- 
tials [41| and the exchange-correlation functionals pa- 
rameterized by Perdew, Burke and Ernzerhof [12] for 
the generalized gradient approximation (GGA). The en- 
ergies were calculated at zero temperature and pressure, 
with spin polarization, and without zero-point motion 
or lattice vibrations. All crystal structures were fully 
relaxed (cell volume and shape and the basis atom coor- 
dinates inside the cell) . Numerical convergence to about 
lmeV/atom was ensured by a high energy cutoff (40% 
higher than the highest energy cutoff for the pseudo- 
potentials of the components) and dense 6000 k-point 
Monkhorst-Pack meshes. For each system, we calculated 
the energies of all the ground state structures reported in 
Refs. [43j and [3] and about 200 additional crystal struc- 
tures. In addition to the 176 structures described in [J] 
these included the prototypes A5, A6, A7, All, B20, Til, 
Thin, LiB-MSl/2 0,111, Au 4 Zr, Ca 7 Ge, NbNi 8 (Pt 8 Ti), 
Ga 2 Hf, WsSis, V4Zn 5 , Ni7Zr 2 , C36 and the complete set 
of hep-superstructures [47J with up to four atoms per cell. 

Some of our results deviate from experimental data 
(43l . |44| . Table U is a summary of the alloy systems ad- 
dressed in this study. On the first column, the alloying 
metals are ordered according to their Mendeleev number 



(or Pettifor's chemical scale) [43, [49|. The reported ex- 
perimental compounds (or lack thereof) are shown in the 
second column and the calculated ones in the third col- 
umn. The compounds are presented with their structure 
Strukturbericht designation or prototype in parentheses 
(unspecified or unknown structures are denoted as Unk.). 
Some of the predicted phases (marked by an asterisk) 
have structures for which no prototype is known and no 
Strukturbericht designation have been given. These new 
prototypes are described in Table [TTJ Appearance of two 
structures in the third column indicates their degeneracy. 
The fourth column gives the calculated energies of the 
structures identified on the convex hull of each binary 
system (3rd column). Energies of reported structures 
(2nd column) that are found to lie above the convex hull 
(missing in the 3rd column) are indicated in parentheses. 
In cases where the reported and calculated structures of 
a compound are different, or two non-degenerate struc- 
tures are reported, the energy difference between them is 
indicated in square parentheses. The calculated convex 
hulls of all the compound forming systems are shown in 
Figs. ©-©. 

In Table U there are 16 systems that are reported in 
the literature as phase separating, i.e. having no com- 
pounds. 13 of these systems are grouped together at the 
top of the table. It is not surprising to see this same 
behavior in a block of the table, because the systems are 
listed by Pettifor's Mendeleev number [49(. The other 
three systems reported to be phase-separating, Hf-Mg, 
Hf-Tl and Hf-Pb, are scattered in the lower rows of the 
table. In these three cases, our calculations show that 
they are in fact compound forming, essentially comple- 
menting the general trend implied by the Pettifor chem- 
ical scale. Stable structures are also found in the binary 
systems of Hf with Ti and Zr (sharing the IVB column 
of the periodic table) and Sc. In five of the six systems 
predicted to be compound forming (excluding Hf-Pb) the 
pure elements have a hep crystal structure and only dis- 
ordered hep solid solutions have been reported for their 
binary alloys with hafnium over the entire range of con- 
centrations. The calculations indicate, however, that in 
four of them stable compounds that are not hep-based 
superstructures should exist at low temperatures. Only 
in the Hf-Sc system, the new compounds are predicted 
to have a hep based structure. Three of these metals 
have very high melting temperatures and the available 
experimental data on their alloys is also limited to high 
temperatures (above 700°C, 800°C, 1000°C for Ti, Zr 
and Sc, respectively [43|, [44j. No data is available for 
Hf-Mg, Hf-Tl and Hf-Pb alloys. 
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TABLE I: Compounds observed in experiments or predicted by ab 
initio calculations in metallic binary alloys of Hf. Structure Struk- 
turbericht designation or prototype are in parentheses (Unk. denotes 
unknown structures). New prototypes are marked by * and described 
Tabic [II] More than one structure may have been reported (2nd 
column) or found with degenerate energies (3rd column). Ai/ are the 
formation enthalpies of the compounds in the present study; parenthe- 
ses denote reported structures (2nd column) that lie above the convex 
hull (missing in the 3rd column). Square parentheses denote energy 
differences between experimentally observed and calculated structures 
or two reported non-degenerate structures. 





Compounds 


AH 




Experiments [43. 441 


Calculations 


mcV/at 


K,Na 






LLBa 






Sr,Ca,Y 






Sc 




Hf 5 Sc* 


-10 






Hf 3 Sc* 


-10 


La 






Zr 


HfZr 2 (C32) 


-22 


Ti 


HfTi 2 (C32) 


-11 


Nb.Ta 






V 


HfV 2 (C15) 


(31) 


Mo 


HfMo 2 (C15) 


HfMo 2 (C15) 


-170 


W 


HfW 2 (C15) 


HfW 2 (C15) 


-171 


Cr 


HfCr 2 (C15, C36) 


HfCr 2 (C15) 


-120[+6] 


Tc 




Hf 3 Tc(Mo 3 Ti* [4]) 


-269 






Hf 2 Tc(C49) 


-357 




HfTc(B2) 


HfTc(B2) 


-482 




HfTc 2 (C14) 


HfTc 2 (C14) 


-362 


Rc 




Hf 3 R.e(Mo 3 Ti* [4]) 


-200 




HfRe (Unk.) [44] 








Hf 2 iRe 25 (Re 25 Zr 2 i) [43] 


Hf 2 iRc 25 (Re 25 Zr 21 ) 


-407 




HfRe 2 (C14) 


HfRe 2 (C14) 


-394 




Hf 5 Re 2 4(Ti 5 Re 2 4) 


Hf 5 Rc 2 4(Ti 5 Re 2 4) 


-252 


Mn 


HfMn 2 (C14, C36) 


HfMn 2 (C14) 


-268[+4] 




Hf 2 Mn(NiTi 2 ) 




(-1091 


Fe 




Fe 5 Hf(C15 b ) 


-178 




Fc 2 Hf(C14, C15) 


Fe 2 Hf(C15) 


-354[+14] 






FeHf(B2) 


-331 




FcHf 2 (NiTi 2 ) 




(-189) 


Os 


Hf 5 40si7(Hf 5 40si 7 ) 




("-SI 9 s ! 




Hf 2 Os(NiTi 2 ) 




(-429) 




HfOs(B2) 


HfOs(B2) 


-707 




HfOs 2 (C14) 




(-402) 


Ru 


HfRu(B2) 


HfRu(B2) 


-819 




HfRu 2 (Unk.) 






Co 


Co 7 Hf 2 (Ni 7 Zr 2 ) 




(-219) 




Co 2 Hf(C15) 


Co 2 Hf(C14) 


-374[+3] 




CoHf(B2) 


CoHf(B33) 


-401[+12] 




CoHf 2 (NiTi 2 ) 


CoHf 2 (C37) 


-314[+23] 


Ir 


Hf 2 Ir(NiTi 2 ) 


Hf 2 Ir(C37) 


-750[+31] 




Hf 5 Ir 3 (D8 8 [43], Ir 5 Zr 3 [44]) 


Hf 5 Ir 3 (Ir 5 Zr 3 ) 


ai 4T^l 41 




Hflr (Unk.) 


Hflr(B27) 


-949 






HfIr 2 (Ga 2 Hf) 


-872 




HfIr 3 (Ll 2 ) 


HfIr 3 (Ll 2 ) 


-800 


Rh 


Hf 2 Rh(NiTi 2 ) 


Hf 2 Rh(CuZr 2 ) 


-632[+12] 




HfRh(B2) 


HfRh(B27) 


-899[+30] 




Hf 3 Rh5(Ge 3 Rh5) 


Hf 3 Rh 5 (Gc 3 Rh5) 


-928 




HfRh 3 (Ll 2 ) 


HfRh 3 (Ll 2 ) 


-762 


Ni 


Hf 2 Ni(C16) 




(-345) 




HfNi(CrB [43], Til [44]) 


HfNi(TlI) 


-542 [+20] 




Hf B Niu(Zr 9 Ptu) 




(-461) 




Hf 7 Ni 10 (Zr 7 Ni 10 ) 




(-520) 




HfNi 3 (BaPb 3 ) 


HfNi 3 (D0 24 ) 


-545 [+4] 




Hf 2 Ni 7 (Ni 7 Zr 2 ) 




(-470) 




HfNi 5 (C15 b ) 




(-350) 


Pt 


Hf 2 Pt(NiTi 2 ) 


Hf 2 Pt(NiTi 2 ) 


-786 




R~fPt(B2, B33 [43], Til [44]) 


HfPt(B33, Til) 


-1153[+165] 





HfPt 3 (Ll 2 , D0 24 ) 


HfPt 3 (D0 24 ) 


-1100[+3] 






HfPts(PtgTi) 


-528 


Pd 


Hf 2 Pd(Cll b [43], CuZr 2 [44]) 


Hf 2 Pd(Cll b , CuZr 2 ) 


-527 




HfPd (Unk.) 


HfPd(B33) 


-682 




Hf 3 Pd 4 (Unk.) 








HfPd 2 (Cll b ) 


HfPd 2 (Cll b ) 


-817 




HfPd 3 (D0 24 , Ll 2 ) 


HfPd 3 (D0 24 ) 


-879[+ll] 






HfPd 5 (HfPd*) 


-635 






HfPd 8 (Pt 8 Ti) 


-430 


Au 


Au 5 Hf [43],Au 4 . 2 Hfo. 8 [44](Dl a ) 




(-407) 




Au 4 Hf(Au 4 Zr,Dl a ) 


Au 4 Hf(Au 4 Zr) 


-414 




Au 3 Hf(D0 a ) 


Au 3 Hf(D0 a ) 


-483 




Au 2 Hf(Cll b ) 


Au 2 Hf(Cll b ) 


-565 






Au 4 Hf 3 (Cu 4 Ti 3 ) 


-563 




AuioHfr (Unk.) 








AuHf(Bll) 


AuHf(Bll) 


-545 




AuHf 2 (Cll b [43], CuZr 2 [44]) 


AuHf 2 (Cll b , CuZr 2 ) 


-440 


Ag 


AgHf(Bll) 


AgHf(Bll) 


-119 




AgHf 2 (Cll b [43], CuZr 2 [44]) 


AgHf 2 (Cll b , CuZr 2 ) 


-122 


Cu 




Cu E Hf(C15i,) 


-129 




Cu 5 lHf 1 4(Ag5lGd 14 ) 




(139) 




Cu 8 Hf 3 (Cu s Hf 3 ) 


Cu 8 Hf 3 (Cu 8 Hf 3 ) 


-173 




Cu 10 Hf 7 (Ni 10 Zr 7 ) 


Cu 10 Hf 7 (Ni 10 Zr 7 ) 


-186 




CuHf 2 (Cll b [43], CuZr 2 [4411 


CuHf 2 (Cll b , CuZr 2 ) 


-166 


Mg 


- 


HfMg(CdTi) 


-7 


Hg 


Hf 2 Hg(Cll b [43], CuZr 2 [44]) 


Hf 2 Hg(Cll b , CuZr 2 ) 


-120 


Cd 


CdHf(Bll [43], CdTi [44]) 


CdHf(CdTi, Bll) 


-88 




CdHf 2 (Cll b [43], CuZr 2 [44]) 


CdHf 2 (/31) 


-87[+18] 


Zn 


Hf 2 Zn(Cll b [43], CuZr 2 [44]) 


Hf 2 Zn(Cll b , CuZr 2 ) 


-175 




HfZn 2 (C36 [43], C15 [44]) 


HfZn 2 (CaIn 2 ) 


-233[+18] 




HfZn 3 (Unk.) 


HfZn 3 (YCd 3 ) 


-204 




HfZn 5 (Unk.) 








HfZn 22 (Zn 22 Zr) 


HfZn 22 (Zn 22 Zr) 


-43 


Be 


Bei 3 Hf(D2 3 ) 


Bei 3 Hf(D2 3 ) 


-173 




Bci 7 Hf 2 (Th 2 Nii 7 ,Th 2 Zm 7 ) 


Bei 7 Hf 2 (Th 2 Zni 7 ) 


-223[+485] 




Be 5 Hf(D2 d ) 


Bc 5 Hf(D2 d ) 


-226 




Bc 2 Hf(C32) 




(-166) 




BcHf(TlI) 




(-75) 


Tl 


Hf 2 Tl(/3 2 ) 


-2 


In 




Hf 3 In(Ll 2 ) 


-145 




Hfln(Llo) 




(-160) 




Hf 3 In 4 (In 4 Ti 3 ) 


Hf 3 In 4 (In 4 Ti 3 ) 


-285 


Al 


Al 3 Hf(D0 22 [43], D0 23 [44]) 


Al 3 Hf(D0 23 ) 


-356[+9] 




Al 2 Hf(C14) 


Al 2 Hf(C14) 


-415 




Al 3 Hf 2 (Al 3 Zr 2 ) 




(-398) 




AlHf(B33 [431. Til [441) 




(-380) 




Al 3 Hf 4 (Al 3 Zr 4 ) 


Al 3 Hf 4 (Al 3 Zr 4 ) 


-374 




Al 2 Hf 3 (Al 2 Zr 3 ) 




(-314) 




Al 3 Hf 5 (D8 s ) 




(-277) 




AlHf 2 (C16) 




(-252) 




AlHf 3 (Unk.) [44] 


AlHf 3 (Ll 2 ) 


-224 


Ga 


Ga 3 Hf(D0 22 ) 




(-382) 




Ga 2 Hf(Ga 2 Hf) 




(-456) 




Ga 3 Hf 2 (Al 3 Zr 2 ) 


Ga 3 Hf 2 (Al 3 Zr 2 ) 


-664 




GaHf(InTh) 




(-519) 




GaioHfn(GeioHon) 




(-495) 




Ga 3 Hf 5 (D8 8 ) 


Ga 3 Hf 5 (D8 8 ) 


-445 




GaHf 2 (C16) 


GaHf 2 (C16) 


-404 


Pb 




Hf 5 Pb* 


-140 






Hf 5 Pb 3 (W 5 Si 3 ) 


-253 


Sn 


Hf 5 Sn 3 (D8 8 ) 


Hf 5 Sn 3 (D8 8 ) 


-380 




Hf 5 Sn 4 (Ga4Ti 5 ) 


Hf5Sn 4 (Ga 4 Ti 5 ) 


-396 




HfSn(B20) 




(-270) 




HfSn 2 (C40) 


HfSn 2 (C40) 


-265 


Bi 


Bi 2 Hf(As 2 Ti) 


Bi 2 Hf(C16) 


-147[+14] 




Bi 9 Hf 8 (V 7 .46Sb 9 ) 




(-166) 




BiHf (Unk.) 


BiHf(Bll) 


-182 






BiHf* 


-169 



4 




5 




FIG. 2: Formation enthalpies of Hf-Mn, Hf-Mo, Hf-Ni, Hf-Os, Hf-Pb, Hf-Pd, Hf-Pt, Hf-Re, Hf-Rh, Hf-Ru, Hf-Sc, Hf-Sn, Hf-Tc, 
Hf-Ti, and Hf-Tl alloys. 
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20 40 60 80 100 

Hf Atomic Percent Zirconium Zr 

FIG. 3: Formation enthalpies of Hf-W, Hf-Zn, and Hf-Zr al- 
loys. 

Following these findings we carried out a more exten- 
sive HT search for additional stable structures in these 
six binary systems. This search included 137 fee struc- 
tures with up to six atoms per cell, 137 bec structures 
with up to six atoms per cell and 333 hep structures 
with up to eight atoms per cell (enumerated in Ref. 
l47t ). This search found lower energy structures in the 
Hf-Sc system but failed to improve on the findings of the 
smaller set search in the other five systems. This demon- 
strates the efficiency and reliability of the HT approach in 
finding phase-stability in binary systems, using a search 
base of about 200 structures per system. More extensive 
searches, on larger sets of derivative structures, seem to 
be useful only in cases where both elements and their 
compounds share a common parent lattice. 

In contrast to the preceding examples, our calculations 
predict phase-separation in the Hf-V system, for which 
the compound HfV2 has been reported at high temper- 
atures (1000-1550 °C) 0,113. Calculations of possible 
structures of this compound show that the lowest ener- 
gies, all positive, are obtained for C14, C36 and C15 in 



that order, whereas the structure reported in the liter- 
ature is C15 1441. This is in agreement with the calcu- 
lation in Ref. [24] that also found C14 as the lowest en- 
ergy structure of these three Laves phases for HfV 2 . The 
compound HfV2 reported in the literature is therefore 
unstable at low temperatures. The discrepancy between 
the computational result of C14 as the minimum energy 
structure (Hf atoms on the sites of a hexagonal diamond 
structure) and the observed high temperature C15 phase 
(Hf on the sites of a diamond structure) may be due to 
vibrational stabilization at high temperature. In fact, 
C15 is only 24 meV/atom above C14, and given that en- 
tropic differences between structures can be of the order 
of 0.1-l.Oke per atom [5(|, very small energy differences 
between the experimentally observed structure and our 
ab initio results could be reversed at elevated tempera- 
ture. The phenomenon is common, e.g. in Ref. l5l| the 
vibrational entropy difference is shown to stabilize the 
fl-A^Cu (C16) phase over the competing A^Cu-fl' phase 
(distortion of (9 C -C1), which has the lowest energy and is, 
therefore, stable at low temperatures. 

In agreement with experimental data, all the other 
transition metals, from columns VIB to IIB of the pe- 
riodic table, have stable compounds with hafnium at low 
temperatures. The situation is especially simple for the 
VIB metals, chromium, molybdenum and tungsten, for 
which the calculations confirm the existence of a single 
compound HfM2 with a C15 structure. Other Hf-Mo 
compounds for which ab initio calculations have been 
previously performed [201 ] are shown to be unstable. The 
Hf-Hg system also exhibits a single compound, Hf 2 Hg, 
where the calculations confirm that the reported struc- 
tures Cllb [43| and CuZr 2 [H| are degenerate. 

The calculations for the Fe-Hf system show the ex- 
istence of a stable compound Fe2Hf with a C15 struc- 
ture, and the C14 structure just slightly less stable (14 
meV/atom). This is consistent with the literature where 
the two structures have been reported [44j at temper- 
atures above 600°C. However, in contrast to the avail- 
able experimental data, the calculation shows that the 
compounds FeHF (B2 structure) and FesHf (C15&) are 
also stable whereas the FeHf2 compound, reported with 
a NiTi2 structure, is unstable at low temperatures. Its 
formation energy is 35 meV/atom higher than the tie 
line FeHf(B2)^Hf(A3) and thus could decompose at 
temperatures lower than 100°C (Fig. [1]). A similar 
behavior is found in the HfMn system where the C14 
structure of HfMn2 is found to be stable, with the C36 
and C15 structures at slightly higher energies (4 and 5 
meV/atom, respectively). The experimentally reported 
Hf2Mn compound (NiTi2 prototype), is again found to 
be unstable at low temperatures, at 30meV/atom above 
the Hf(A3)<^HfMn 2 (C14) tie line. 

In the Hf-Ru and Hf-Os systems, a stable B2 struc- 
ture of HfRu and HfOs is found in agreement with 
the experiments (44|. The compounds Hf20s, Hfs40si7 
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and HfOs2 reported from experimental data are found 
to be unstable. No stable structure was found for 
the HfRu2 compound, whose existence was suspected 
in some experiments [Hj]. In the Hf-Re system rcf. 
[H| reports a Hf 2 i Re 2 s compound (prototype Re25Zr 2 i) 
whereas ref. [43| reports a HfRe compound of unspec- 
ified prototype. Our calculations confirm the existence 
of the Hf2iRe 2 5 compound but find no stable structure 
for HfRh. The lowest HfRh structure, B2, lies approxi- 
mately 50meV/atom above the convex hull (Table U and 
Fig. [2]). For the Co-Hf system, our calculations confirm 
the existence of the stable compounds Co2Hf, CoHf and 
CoHf 2 at low temperatures, albeit with structures differ- 
ent from those reported in experiments. The reported 
NiyZr2 structure of Co7Hf 2 lies above the convex hull 
(Table Hand Fig. [TJ. 

The Cu-Hf system is the only transition metal-hafnium 
system previously studied by ab initio methods, also us- 
ing the VASP software but with different potentials than 
in our study (Vanderbilt-type ultrasoft pseudopotentials) 
[37l | . 28 different compound structures have been calcu- 
lated, of which four, CusHf, CusHf 3 , CuioHf? and CuHf 2 
were identified as stable. Our calculations also indicate 
that these four compounds are stable. There is no exper- 
imental information about the sructure of CusHf [44j. 
We obtain the AuBes structure in agreement with Rcf. 
(37j . The CuZr2 structure, reported in experiments for 
CuHf2 [44], and not calculated in ref. [33] is degener- 
ate with the reported Cllf, structure. The ground state 
convex hull is asymmetric and skewed towards the Cu 
side, in agreement with the general conclusions of Rcf. 
[33] (Fig. [[]). The Al-Hf system was studied using the 
same methodology 35]. The convex hull in this case 
is defined by four compounds Al 3 Hf(D0 2 3), Al 2 Hf(C14), 
Al 3 Hf4(Al 3 Zr4) and AlHf 3 (Ll 2 ), with Al 3 Hf 2 (Al 3 Zr 2 ) 
and AlHf(B33) slightly above the convex hull. Our cal- 
culations reproduce these results and show in addition 
that the AlHf structure Til, reported in the experimen- 
tal data but not studied in Ref. 135], is degenerate with 
B33. Other compounds reported for this system [3] are 
found to be significantly above the convex hull and should 
therefore be unstable at low temperatures (see Table Q] 
and Fig.CQ). 

Our HT study uncovers a few compounds on which no 
data is available in the experimental literature [HI, [44| in 
the binary systems of Hf with Au, Bi, In, Pd, Pt, Re, and 
Tc. A few reported compounds are found to lie above the 
convex hull of their respective systems (Hf with Be, Ga, 
Ni and Sn) and are thus predicted to be unstable at low 
temperatures. (Table HJ 

In most cases where two structures are reported in 
the experimental literature, the calculations show that 
they are indeed degenerate or the energy differences are 
small, 30 meV/atom or less (square parentheses in Ta- 
ble |T|. In only two cases these structures differ con- 
siderably. The B2 structure of HfPt is 165meV/atom 



higher than the degenerate B33/T1I groundstate. No 
phase diagram is available for the Hf-Pt system, but 
the experimenatal data indicates that the B2 structure 
might have been observed in a non-stoichiomctric mix- 
ture (43| . The calculations seem to provide an indirect 
confirmation of this observation, showing that B2 is not 
a stable structure of the stoichiometric mixture. The 
TI12N117 structure of Bei7Hf2 is 485meV/atom higher 
than the Th ^Zn -i? g roundstate (sometimes also denoted 
as Bei7Nb2 [43L El]). As discussed in Ref. 0, it has 
not been determined whether Bei7Hf2 consists of both 
forms or just one of them. The large energy difference 
obtained in our calculations indicates that one structure 
is much more likely than the other and the compound 
should consist exclusively of the Tti2Zri7 prototype. 

In conclusion, a systematic and comprehensive ab 
initio study of phase stability is carried out for the 
hafnium intermetallic systems. The total energies of 
about 200 intermetallic compounds have been calculated 
for 44 hafnium-metal systems and some interesting de- 
viations from published experimental data were found. 
In particular, the calculations predict the existence of 
stable compounds in six hafnium intermetallic systems 
previously believed to be phase-separating. On the Pet- 
tifor chemical scale, three of these are isolated phase- 
separating systems in a cluster of compound-forming 
ones, and the discrepency is likely due to lack of exper- 
imental data. Hence, these predictions nicely comple- 
ment the trend indicated by the empirical scale. A few 
new compounds are predicted in binary systems of Hf 
and other metals, and some compounds reported in the 
literature are shown to be unstable at low temperatures. 

A detailed understanding of Hf alloys is crucial for a 
better realization of its potential as an alloying agent in 
currently available applications and in developing new 
ones. The picture of Hf alloys that emerges from this 
study is quite different from that depicted by current ex- 
perimental data. It should be emphasized that we con- 
sider the alloys to be in thermodynamical equilibrium, 
which can be difficult to reach at low temperatures due 
to slow kinetics. At higher temperatures, configurational 
disorder and vibrational entropic promotion might desta- 
bilize the predicted compounds. The theoretical predic- 
tions presented here should therefore motivate research 
for their experimental validation and provide useful guid- 
ance to future studies of these alloys. 
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TABLE II: Geometry of new prototypes in our study. Positions are given as unrelaxed positions in the parent lattice. 



System 


Hf 3 Sc* 


Hf 5 Sc* 


BiHf* 


Mo 3 Ti* 


Hf 5 Pb* 


HfPd* 


Superlatticc 


HCP A 6 B 2 


HCP A 5 B! 


HCP A 2 B 4 


BCC AB 3 
[|] 


FCC AsBj 


none 


Lattice 


Orthorhombic 


Hexagonal 


Monoclinic 


Orthorhombic 


Tetragonal 


Orthorhombic 


Space Group 


Cmcm #63 


P62m #189 


C2/m # 12 


Immra #71 


P4/mmm #123 


Cmmm #65 


Pearson symbol 


oC16 


hP6 


mC12 


oI8 


tP6 


oC12 


Primitive 














vectors (cart.) 














ai/a 


(2,0,0) 


(3/2, %/3/2,0) 


(-1,0,0) 


(3/2, 1/2, -1/2) 


(1/2,1/2,0) 


(-2.01059,5.93236,0) 


a 2 /a 


(1, V3, 0) 


(3/2,3\/3/2,0) 


(1/2, V3/2, \/8/3) 


(1/2,3/2,1/2) 


(0,3,3) 


(-2.01059, -5.93236,0) 


a 3 /a 


(0,0, s/8/3) 


(0,0, vWa) 


(o, V3. -vWs) 


(-1/2, -3/2, 1/2) 


(1/2,5/2,3) 


(0,0,4.02718) 


Atomic 














positions (fract.) 














Al 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(0.02752, 0.97247, 0) 


A2 


(0, 1/2,0) 


(0, 1/3,0) 


(7/18, 7/95/18) 




(0, 1/6,0) 




A3 


(1/2,0,0) 


(0,2/3,0) 






(0, 1/3,0) 




A4 


(1/6, 1/6, 1/2) 


(1/3, 0, 1/2) 






(0, 1/2,0) 




A5 


(1/6, 2/3. 1/2) 


(1/3, 1/3. 1/2) 






(0, 2/3, 0) 




A6 


(2/3, 1/6, 1/2) 












Bl 


(1/2, 1/2,0) 


(1/3,2/3, 1/2) 


(1/3,2/3,2/3) 


(1/4,3/4,1/2) 


(0,5/6,0) 


(0.19430,0.80569,0.5) 


B2 


(2/3,2/3,1/2) 




(2/3, 1/3, 1/3) 


(1/2,1/2,0) 




(0.36452,0.63547,0) 


B3 






(13/18,4/9, 17/18) 


(3/4,1/4,1/2) 




(0.52795,0.47204,0.5) 


B4 






(1/18, 1/9, 11/18) 






(0.69061,0.30938,0) 


B5 












(0.86173,0.13826,0.5) 


AFLOW label 


"471" 


"473" 


"475" 


"81" 


"477" 


"479" 


[39] 
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